Tarea 3: Clasificación de texto con RNN

CC6204 - Deep Learning


Cuerpo Docente:¶

  • Profesor: Iván Sipiran
  • Ayudantes:
    • Camila Figueroa Acevedo
    • Gustavo Santelices
    • Sofia Capibara Chávez Bastidas
    • Victor Faraggi V. ### Estudiante:
  • Maximiliano Varas González (maximilianovarasg@gmail.com)

En esta tarea van a crear una red neuronal que clasifique mensajes como spam o no spam. Lo primero es descargar la data:

In [1]:
!mkdir data
!python -m wget https://www.ivan-sipiran.com/downloads/spam.csv
!move *.csv data/
Ya existe el subdirectorio o el archivo data.
Saved under spam.csv
C:\Users\mvarasg\Documents\Universidad\CC6204 - Deep Learning\Tareas\T3\spam.csv
Se han movido         1 archivos.

Los datos vienen en un archivo CSV que contiene dos columnas "text" y "label". La columna "text" contiene el texto del mensaje y la columna "label" contiene las etiquetas "ham" y "spam". Un mensaje "ham" es un mensaje que no se considera spam.

Tarea¶

El objetivo de la tarea es crear una red neuronal que clasifique los datos entregados. Para lograr esto debes:

  • Implementar el pre-procesamiento de los datos que creas necesario.
  • Particionar los datos en 70% entrenamiento, 10% validación y 20% test.
  • Usa los datos de entrenamiento y valiadación para tus experimentos y sólo usa el conjunto de test para reportar el resultado final.

Para el diseño de la red neuronal puedes usar una red neuronal recurrente o una red basada en transformers. El objetivo de la tarea no es obtener el performance ultra máximo, sino entender qué decisiones de diseño afectan la solución de un problema como este. Lo que si es necesario (como siempre) es que discutas los resultados y decisiones realizadas.


Pre-procesamiento de data¶

Importando librerías:

In [2]:
import pandas as pd
import numpy as np
import random
import torch
from torch.utils.data import TensorDataset, DataLoader
import matplotlib.pyplot as plt
import plotly.express as px
import torch.nn as nn
from sklearn.metrics import confusion_matrix
from sklearn.model_selection import train_test_split
import seaborn as sns
from string import punctuation
import plotly
plotly.offline.init_notebook_mode()
import os 
# from sklearn.utils.multiclass import unique_labels
# from sklearn.metrics import classification_report


import time
SEED = 1234
random.seed(SEED)
np.random.seed(SEED)
torch.manual_seed(SEED)
torch.cuda.manual_seed(SEED)
torch.backends.cudnn.deterministic = True

Se carga la data en un DataFrame:

In [3]:
data = pd.read_csv('./data/spam.csv')
print(f'Las dimensiones del Dataframe son: {data.shape}.')
Las dimensiones del Dataframe son: (5572, 2).

Veamos si hay elementos nulos en ambas columnas

In [4]:
data.isna().sum()
Out[4]:
text     1
label    1
dtype: int64

En efecto, hay al menos una entrada en la columna text y en label.

In [5]:
data[data.isna().any(axis=1)]
Out[5]:
text label
2115 Well I wasn't available as I washob nobbing wi... NaN
3035 NaN -) ok. I feel like john lennon.

Hay elementos NaN en filas diferentes, por lo que se eliminan estas dos.

In [6]:
data = data.dropna()

Veamos que elementos hay en la columna label, debería solo estar "spam" y "ham".

In [7]:
# print(data['label'].unique())
print(f"En total hay {len(data['label'].unique())} distintas labels, cuando deberían ser dos.")
print(f"Las entradas distintas de 'spam' y 'ham' ({len(data[(data['label']!='spam') & (data['label']!='ham')])}) son: ")

data[(data['label']!='spam') & (data['label']!='ham')]
En total hay 43 distintas labels, cuando deberían ser dos.
Las entradas distintas de 'spam' y 'ham' (207) son: 
Out[7]:
text label
44 Great! I hope you like your man well endowed. ... #&gt
55 Do you know what Mallika Sherawat did yesterda... URL&gt
78 Does not operate after &lt #&gt
200 I sent you &lt #&gt
202 Your account has been refilled successfully by... DECIMAL&gt
... ... ...
5500 Love has one law Make happy the person you love. In the same w...
5504 Wait . I will msg after &lt #&gt
5513 Yes. Please leave at &lt #&gt
5557 No. I meant the calculation is the same. That ... #&gt
5559 if you aren't here in the next &lt #&gt

207 rows × 2 columns

Se procede a eliminar los elementos que no tiene bien la etiqueta.

In [8]:
data = data[(data['label'] == 'spam') | (data['label'] == 'ham')]

Se procede a eliminar los signos de puntuación, entradas vacías y '\n' de las entradas. Luego, se pasa todo el texto a minúsculas.

In [9]:
print(punctuation)
data['text'] = data['text'].str.lower()
data['text'] = data['text'].str.replace(f'[{punctuation}]', '', regex=True)
data = data[data['text']!=' '] #Se eliminan elementos con strings vacíos.
data['text'] = data['text'].str.replace('\n', ' ')
print(f'Al eliminar elementos NaN y vacíos la nueva dimensión de data es: {data.shape}.')
!"#$%&'()*+,-./:;<=>?@[\]^_`{|}~
Al eliminar elementos NaN y vacíos la nueva dimensión de data es: (5361, 2).

Se pasan todas las palabras de la columna text a una lista y se hace un split de esta para tener todas las palabras:

In [10]:
lista_words = data['text'].tolist()
lista_words = [word for instancia in lista_words for word in instancia.split()]
print(f'En total se tienen {len(lista_words)} palabras.')
En total se tienen 81325 palabras.

Se crea el vocabulario para utilizar más adelante con el Embedding.

In [11]:
from collections import Counter

counts = Counter(lista_words) #Diccionario donde la llave es la palabra y el valor es la frecuencia.
vocab = sorted(counts, key=counts.get, reverse=True) #Lista de palabras, ordenadas de mayor a menor.
vocab_to_int = {word: ii for ii, word in enumerate(vocab, 1)} #Diccionario con palabra como key y la posición de las más frecuentes como value.
print('Palabras únicas:', len(vocab_to_int))
Palabras únicas: 9328

Se agrega una nueva columna con un representación de lista de números del textp. Se utiliza el vocabulario.

In [12]:
data['text_num'] = data['text'].apply(lambda x: [vocab_to_int.get(word, word) for word in x.split()]) #Nueva columna con la lista con la representación numerica de cada palabras.
data.head(5)
Out[12]:
text label text_num
0 go until jurong point crazy available only in ... ham [46, 436, 4264, 773, 728, 729, 64, 8, 1194, 88...
1 ok lar joking wif u oni ham [48, 301, 1356, 418, 6, 1782]
2 free entry in 2 a wkly comp to win fa cup fina... spam [45, 437, 8, 22, 4, 730, 871, 1, 175, 1783, 11...
3 u dun say so early hor u c already then say ham [6, 217, 143, 24, 353, 2828, 6, 156, 140, 58, ...
4 nah i dont think he goes to usf he lives aroun... ham [921, 2, 50, 99, 70, 438, 1, 922, 70, 1786, 21...

Sigue reemplazar la columna label por valores numéricos. Se considera que spam:1 y ham:0.

In [13]:
data['label_encoded'] = data['label'].replace({'spam':1, 'ham':0})

Padding¶

Consiste en agregar o eliminar palabras con el fin de que la entrada a la red sea siempre la misma. Para decidir el tamaño de las secuencias de texto, se obtiene información estadística del largo de texto del dataset.

In [14]:
from statistics import median, mode
data['largo_text'] = data['text_num'].apply(len)
print("Largo promedio de las listas:", data['largo_text'].mean())
print("Mediana del largo de las listas:", median(data['largo_text']))
print("Moda del largo de las listas:", mode(data['largo_text']))
print('Max:', max(data['largo_text']))
print('Min:', min(data['largo_text']))




fig = px.histogram(data, x='largo_text', nbins=200, title='Distribución largo texto',marginal='box')

fig.show()
Largo promedio de las listas: 15.16974445066219
Mediana del largo de las listas: 12
Moda del largo de las listas: 6
Max: 171
Min: 1

En base a esto se observa que la gran mayoría de los textos tiene largo bajo, específicamente largo 6. Sin embargo, se estima que definir un largo de secuencia tan pequeño puede ser muy perjudicial, dado que gran cantidad de textos tienen más información que 6 palabras. Se considera aumentar el largo de secuencia no representa un salto sustancial en el costo de entrenamiento, siempre y cuando este valor no sea muy grande. Viendo la distribución, se decide tomar el largo de secuencia igual al upper fence (44), dado que engloba a la mayoría de datos y porque cadenas de textos más grandes son outliers o poco comunes.

Se define la siguiente estrategia:

  • Se define un tamaño para las secuencias de texto: largo_seq
  • Si el largo de una secuencia es mayor a largo_seq, se corta y se conservan las primeras largo_seq palabras.
  • Si el largo de una secuencia es menor a largo_seq, se rellena con 0's al principio y se dejan al final las palabras de la secuencia, con el fin de tener una secuencia de largo largo_seq.

Se realiza lo anterior y se agrega la lista con padding en una nueva columna del DataFrame.

In [15]:
def padding_seq(data, largo_seq):
    if isinstance(data, list):
        features = np.zeros((1, largo_seq), dtype=int)
        features[0, -len(data):] = np.array(data)[:largo_seq]
        return features
    elif isinstance(data, pd.DataFrame):
        df = data.copy()
        def padd_row(row, largo_seq):
            if row['largo_text'] > largo_seq: # Se trunca
                return row['text_num'][:largo_seq]
            elif row['largo_text'] < largo_seq: # Se rellena con 0's
                return [0] * (largo_seq - row['largo_text']) + row['text_num']
            else: # Largo igual
                return row['text_num']
        df['text_num_padded'] = df.apply(padd_row, args=(largo_seq,), axis=1)
        return df
    else:
        raise ValueError("El argumento debe ser una lista de números o un DataFrame.")

data_padding = padding_seq(data, 44)
data_padding
Out[15]:
text label text_num label_encoded largo_text text_num_padded
0 go until jurong point crazy available only in ... ham [46, 436, 4264, 773, 728, 729, 64, 8, 1194, 88... 0 20 [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...
1 ok lar joking wif u oni ham [48, 301, 1356, 418, 6, 1782] 0 6 [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...
2 free entry in 2 a wkly comp to win fa cup fina... spam [45, 437, 8, 22, 4, 730, 871, 1, 175, 1783, 11... 1 28 [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...
3 u dun say so early hor u c already then say ham [6, 217, 143, 24, 353, 2828, 6, 156, 140, 58, ... 0 11 [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...
4 nah i dont think he goes to usf he lives aroun... ham [921, 2, 50, 99, 70, 438, 1, 922, 70, 1786, 21... 0 13 [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...
... ... ... ... ... ... ...
5567 this is the 2nd time we have tried 2 contact u... spam [40, 9, 5, 384, 63, 38, 17, 526, 22, 185, 6, 6... 1 30 [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 40,...
5568 will ì b going to esplanade fr home ham [35, 106, 192, 73, 1, 1881, 843, 80] 0 8 [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...
5569 pity was in mood for that soany other suggest... ham [9325, 59, 8, 1191, 12, 19, 9326, 224, 9327] 0 9 [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...
5570 the guy did some bitching but i acted like id ... ham [5, 521, 108, 109, 9328, 26, 2, 4206, 57, 391,... 0 26 [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...
5571 rofl its true to its name ham [2509, 42, 563, 1, 42, 260] 0 6 [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...

5361 rows × 6 columns


Particionamiento de datos¶

Se divide en 70% para el conjunto de train, 15% para val y 15% para test.

In [16]:
proporcion_train = 0.7
proporcion_val = 0.15
proporcion_test = 0.15

# Dividimos el DataFrame en train, val y test
data_train, data_temp = train_test_split(data_padding, test_size=1 - proporcion_train, random_state=SEED)
data_val, data_test = train_test_split(data_temp, test_size=proporcion_test / (proporcion_val + proporcion_test), random_state=SEED)

# Imprime la longitud de los conjuntos para verificar la división
print("Número de ejemplos en train:", len(data_train))
print("Número de ejemplos en val:", len(data_val))
print("Número de ejemplos en test:", len(data_test))
Número de ejemplos en train: 3752
Número de ejemplos en val: 804
Número de ejemplos en test: 805
In [17]:
print("\t\t\tFeatures:")
print("Train set: \t\t{}".format(np.array(data_train['text_num_padded'].to_list()).shape),
      "\nValidation set: \t{}".format(np.array(data_val['text_num_padded'].to_list()).shape),
      "\nTest set: \t\t{}".format(np.array(data_test['text_num_padded'].to_list()).shape))
			Features:
Train set: 		(3752, 44) 
Validation set: 	(804, 44) 
Test set: 		(805, 44)

Se crean los dataloaders:

In [18]:
#Se pasa a tensores
train_x = torch.tensor(np.array(data_train['text_num_padded'].to_list()), dtype=torch.long)  # Usar dtype=torch.long
train_y = torch.tensor(data_train['label_encoded'].to_list(), dtype=torch.float32)

val_x = torch.tensor(np.array(data_val['text_num_padded'].to_list()), dtype=torch.long)  # Usar dtype=torch.long
val_y = torch.tensor(data_val['label_encoded'].values, dtype=torch.float32)

test_x = torch.tensor(np.array(data_test['text_num_padded'].to_list()), dtype=torch.long)  # Usar dtype=torch.long
test_y = torch.tensor(data_test['label_encoded'].values, dtype=torch.float32)

#TensorDataset
train_data = TensorDataset(train_x, train_y)
valid_data = TensorDataset(val_x, val_y)
test_data = TensorDataset(test_x, test_y)

batch_size = 40

#DataLoaders
train_loader = DataLoader(train_data, shuffle=True, batch_size=batch_size, drop_last=True, generator=torch.Generator().manual_seed(SEED))
valid_loader = DataLoader(valid_data, shuffle=True, batch_size=batch_size, drop_last=True, generator=torch.Generator().manual_seed(SEED))
test_loader = DataLoader(test_data, shuffle=False, batch_size=batch_size, drop_last=True, generator=torch.Generator().manual_seed(SEED))

Se revisa si se entrena en CPU o GPU:

In [19]:
# Chequear si tenemos GPU
train_on_gpu=torch.cuda.is_available()

if(train_on_gpu):
    print('Training on GPU.')
else:
    print('No GPU available, training on CPU.')
Training on GPU.

Se crea la arquitectura de la RNN utilizando como base la red del laboratorio 11. A grandes rasgos, consiste en un Embedding, un LSTM, dropout, capa fc (lineal) y una sigmoide.

In [20]:
#Creamos la red neuronal

class SpamRNN(nn.Module):
    def __init__(self, vocab_size, output_size, embedding_dim, hidden_dim, n_layers, drop_prob=0.5):
        super(SpamRNN, self).__init__()

        self.output_size = output_size
        self.n_layers = n_layers
        self.hidden_dim = hidden_dim
        
        # Capas embedding y LSTM
        self.embedding = nn.Embedding(vocab_size, embedding_dim)
        self.lstm = nn.LSTM(embedding_dim, hidden_dim, n_layers,
                            dropout=drop_prob, batch_first=True)
        
        # dropout
        self.dropout = nn.Dropout(drop_prob)
        
        # Capa lineal y sigmoide
        self.fc = nn.Linear(hidden_dim, output_size)
        self.sig = nn.Sigmoid()

    def forward(self, x, hidden):
        
        embeds = self.embedding(x)
        lstm_out, hidden = self.lstm(embeds, hidden)
                
        #Tomamos solo el último valor de salida del LSTM
        lstm_out = lstm_out[:,-1,:]
                
        # dropout y fully-connected
        out = self.dropout(lstm_out)
        out = self.fc(out)
               
        # sigmoide
        sig_out = self.sig(out)
                  
        # retornar sigmoide y último estado oculto
        return sig_out, hidden
    
    
    def init_hidden(self, batch_size):
        # Crea dos nuevos tensores con tamaño n_layers x batch_size x hidden_dim,
        # inicializados a cero, para estado oculto y memoria de LSTM
        weight = next(self.parameters()).data
        
        if(train_on_gpu):
          hidden = (weight.new(self.n_layers, batch_size, self.hidden_dim).zero_().cuda(),
                   weight.new(self.n_layers, batch_size, self.hidden_dim).zero_().cuda())
        else:
          hidden = (weight.new(self.n_layers, batch_size, self.hidden_dim).zero_(),
                   weight.new(self.n_layers, batch_size, self.hidden_dim).zero_())
        # print("Hidden size:", hidden[0].size()) 
        return hidden

Función para calcular tiempo transcurrido:

In [21]:
def time_elapsed(ti, tf):
    t = tf - ti
    h, m = divmod(t, 3600)
    m, s = divmod(m, 60)
    time_str = "Time Elapsed:"
    if h > 0:
        time_str += f" {int(h)} hrs" if h > 1 else f" {int(h)} hr"
    if m > 0:
        time_str += f" {int(m)} mins" if m > 1 else f" {int(m)} min"
    if s > 0:
        time_str += f" {round(s, 2)} segs" if s != 1 else f" {round(s, 2)} seg"
    print(time_str)
    return h, m, s

Se define la función para entrenar la red.

In [22]:
def train(red, epochs, lr=0.001, batch_size=44, clip=5):
    ti = time.time()
    criterion = nn.BCELoss()
    optimizer = torch.optim.Adam(red.parameters(), lr=lr)

    # Listas para almacenar las métricas en cada época
    train_losses = []
    train_accs = []
    val_losses = []
    val_accs = []

    # Enviar red al GPU
    if train_on_gpu:
        red.cuda()

    red.train()

    for e in range(epochs):
        # Inicializar estado oculto
        h = red.init_hidden(batch_size)

        # Variables para calcular las métricas en cada época
        epoch_train_losses = []
        epoch_train_accs = []
        epoch_val_losses = []
        epoch_val_accs = []

        # Bucle para batchs
        for inputs, labels in train_loader:
            if train_on_gpu:
                inputs, labels = inputs.cuda(), labels.cuda()

            # Crear nuevas variables para estados ocultos
            h = tuple([each.data for each in h])

            red.zero_grad()

            # Hacer pasada forward
            output, h = red(inputs, h)
            loss = criterion(output.squeeze(), labels.float())
            loss.backward()
            # Gradient clipping
            nn.utils.clip_grad_norm_(red.parameters(), clip)
            optimizer.step()

            # Calcular precisión en entrenamiento
            pred = torch.round(output.squeeze())
            correct_tensor = pred.eq(labels.float().view_as(pred))
            accuracy = correct_tensor.sum().item() / correct_tensor.numel()
            epoch_train_accs.append(accuracy)
            epoch_train_losses.append(loss.item())

        # Calcular promedio de pérdida y precisión en entrenamiento
        mean_train_loss = np.mean(epoch_train_losses)
        mean_train_acc = np.mean(epoch_train_accs)
        train_losses.append(mean_train_loss)
        train_accs.append(mean_train_acc)

        # Validation loss
        val_h = red.init_hidden(batch_size)
        red.eval()
        for inputs, labels in valid_loader:
            val_h = tuple([each.data for each in val_h])
            if train_on_gpu:
                inputs, labels = inputs.cuda(), labels.cuda()

            output, val_h = red(inputs, val_h)
            loss = criterion(output.squeeze(), labels.float())
            epoch_val_losses.append(loss.item())

            # Calcular precisión en validación
            pred = torch.round(output.squeeze())
            correct_tensor = pred.eq(labels.float().view_as(pred))
            accuracy = correct_tensor.sum().item() / correct_tensor.numel()
            epoch_val_accs.append(accuracy)

        # Calcular promedio de pérdida y precisión en validación
        mean_val_loss = np.mean(epoch_val_losses)
        mean_val_acc = np.mean(epoch_val_accs)
        val_losses.append(mean_val_loss)
        val_accs.append(mean_val_acc)

        red.train()

    # Imprimir métricas de la época
        print(f"Época: {e + 1}/{ epochs}")
        print(f"Train Loss: {round(mean_train_loss,6)} \t Train Accuracy: {round(mean_train_acc,3)}")
        print(f"Val Loss: {round(mean_val_loss,6)} \t Val Accuracy: {round(mean_val_acc,3)}")
        print(7*"--------")
    tf = time.time()
    h, m, s = time_elapsed(ti, tf)
    return {
        'red_entrenada': red,
        'train_losses': train_losses,
        'train_accs': train_accs,
        'val_losses': val_losses,
        'val_accs': val_accs,
        'n_epochs': epochs,
    }

Se define la red dada la arquitectura y se entrena.

In [23]:
vocab_size = len(vocab_to_int) + 1 # +1 for zero padding + our word tokens
output_size = 1
embedding_dim = 100 #Se define un embedding pequeño dado que son pocas palabras en el vocabulario.
hidden_dim = 256
n_layers = 2
net = SpamRNN(vocab_size, output_size, embedding_dim, hidden_dim, n_layers)
dict_red1 = train(net, epochs = 20, lr = 0.001, batch_size = batch_size, clip = 5)
Época: 1/20
Train Loss: 0.213428 	 Train Accuracy: 0.929
Val Loss: 0.085864 	 Val Accuracy: 0.981
--------------------------------------------------------
Época: 2/20
Train Loss: 0.261135 	 Train Accuracy: 0.905
Val Loss: 0.361502 	 Val Accuracy: 0.857
--------------------------------------------------------
Época: 3/20
Train Loss: 0.133433 	 Train Accuracy: 0.955
Val Loss: 0.094756 	 Val Accuracy: 0.975
--------------------------------------------------------
Época: 4/20
Train Loss: 0.05997 	 Train Accuracy: 0.985
Val Loss: 0.054669 	 Val Accuracy: 0.986
--------------------------------------------------------
Época: 5/20
Train Loss: 0.046666 	 Train Accuracy: 0.987
Val Loss: 0.052205 	 Val Accuracy: 0.983
--------------------------------------------------------
Época: 6/20
Train Loss: 0.043256 	 Train Accuracy: 0.987
Val Loss: 0.066617 	 Val Accuracy: 0.985
--------------------------------------------------------
Época: 7/20
Train Loss: 0.024785 	 Train Accuracy: 0.995
Val Loss: 0.057388 	 Val Accuracy: 0.986
--------------------------------------------------------
Época: 8/20
Train Loss: 0.013696 	 Train Accuracy: 0.997
Val Loss: 0.070463 	 Val Accuracy: 0.985
--------------------------------------------------------
Época: 9/20
Train Loss: 0.009281 	 Train Accuracy: 0.997
Val Loss: 0.075288 	 Val Accuracy: 0.984
--------------------------------------------------------
Época: 10/20
Train Loss: 0.003934 	 Train Accuracy: 0.999
Val Loss: 0.094247 	 Val Accuracy: 0.982
--------------------------------------------------------
Época: 11/20
Train Loss: 0.001502 	 Train Accuracy: 1.0
Val Loss: 0.080348 	 Val Accuracy: 0.986
--------------------------------------------------------
Época: 12/20
Train Loss: 0.000567 	 Train Accuracy: 1.0
Val Loss: 0.095039 	 Val Accuracy: 0.985
--------------------------------------------------------
Época: 13/20
Train Loss: 0.000314 	 Train Accuracy: 1.0
Val Loss: 0.08988 	 Val Accuracy: 0.986
--------------------------------------------------------
Época: 14/20
Train Loss: 0.000175 	 Train Accuracy: 1.0
Val Loss: 0.092496 	 Val Accuracy: 0.986
--------------------------------------------------------
Época: 15/20
Train Loss: 0.000135 	 Train Accuracy: 1.0
Val Loss: 0.095437 	 Val Accuracy: 0.985
--------------------------------------------------------
Época: 16/20
Train Loss: 8.9e-05 	 Train Accuracy: 1.0
Val Loss: 0.098536 	 Val Accuracy: 0.985
--------------------------------------------------------
Época: 17/20
Train Loss: 7.6e-05 	 Train Accuracy: 1.0
Val Loss: 0.100268 	 Val Accuracy: 0.984
--------------------------------------------------------
Época: 18/20
Train Loss: 7.3e-05 	 Train Accuracy: 1.0
Val Loss: 0.100868 	 Val Accuracy: 0.984
--------------------------------------------------------
Época: 19/20
Train Loss: 5.4e-05 	 Train Accuracy: 1.0
Val Loss: 0.103054 	 Val Accuracy: 0.984
--------------------------------------------------------
Época: 20/20
Train Loss: 4.6e-05 	 Train Accuracy: 1.0
Val Loss: 0.104745 	 Val Accuracy: 0.984
--------------------------------------------------------
Time Elapsed: 28.57 segs

Función para graficar el Loss y el Accuracy en el conjunto de entrenamiento y validación en las épocas, resultante del entrenamiento.

In [24]:
def graficar(dict_model):
    train_losses = dict_model['train_losses']
    val_losses = dict_model['val_losses']
    train_accs = dict_model['train_accs']
    val_accs = dict_model['val_accs']
    n_epoch = dict_model['n_epochs']
    epoch = range(1,len(train_losses)+1)

    fig, axes = plt.subplots(nrows=1, ncols=2, figsize=(12, 5))  # Dos gráficos en una fila
    axes[0].plot(epoch, train_losses, label='Train Losses', color='b')
    axes[0].plot(epoch, val_losses, label='Val Losses', color='r')
    axes[0].set_title('Función de Pérdidas')
    axes[0].set_xlabel('Épocas')
    axes[0].set_ylabel('Loss')
    axes[0].legend()
    axes[1].plot(epoch, train_accs, label='Train_Acc', color='b')
    axes[1].plot(epoch, val_accs, label='Val_Acc', color='r')


    axes[1].set_title('Accuracy')
    axes[1].set_xlabel('Épocas')
    axes[1].set_ylabel('Accuracy')
    axes[1].legend()
    plt.suptitle(f'Entrenamiento con {n_epoch} epocas.')
    plt.tight_layout()
    plt.show()
In [25]:
graficar(dict_red1)

Los graficos muestran como cambia el Loss y el accuracy en train y val. Se puede ver como al principio la función de loss aumenta y empeora el accuracy, pero se soluciona más adelante. Notar que la diferencia entre el accuracü de train y val es menor al 5%, por lo que se considera que hay poco sobreajuste. Sin embargo, se puede ver como la red comienza a hacerlo, basta ver el loss tiende a cero en train y en val comienza a subir.

Guardar modelo:

In [35]:
if os.path.exists('./SpamRNN.pt'):
    # Si el archivo existe, cargar el modelo
    print('Ya existe el modelo en la ruta.')
    #model = torch.load('./SpamRNN.pt')
else: #guardar
    torch.save(dict_red1['red_entrenada'], './SpamRNN.pt')

Función para evaluar la red en el conjunto de prueba y grafica una matriz de confusión:

In [27]:
def evaluar_red(dict_red, test_loader=test_loader):
    ti = time.time()

    red = dict_red['red_entrenada']
    criterion = nn.BCELoss()

    test_losses = []  # Track loss
    num_correct = 0
    all_preds = []
    all_labels = []

    # Inicializar estado oculto
    h = red.init_hidden(batch_size)  # Asegúrate de tener la variable batch_size definida

    red.eval()
    for inputs, labels in test_loader:

        h = tuple([each.data for each in h])

        if train_on_gpu:
            inputs, labels = inputs.cuda(), labels.cuda()

        output, h = red(inputs, h)

        test_loss = criterion(output.squeeze(), labels.float())
        test_losses.append(test_loss.item())

        # Convertir probabilidades a clases (0,1)
        pred = torch.round(output.squeeze())

        # Comparar predicciones a labels
        correct_tensor = pred.eq(labels.float().view_as(pred))
        correct = np.squeeze(correct_tensor.cpu().numpy()) if train_on_gpu else np.squeeze(correct_tensor.numpy())
        num_correct += np.sum(correct)

        all_preds.extend(pred.detach().cpu().numpy())
        all_labels.extend(labels.cpu().numpy())

    # Calcular estadísticas
    test_loss_avg = np.mean(test_losses)
    test_accuracy = num_correct / len(test_loader.dataset)

    # Matriz de confusión
    cm = confusion_matrix(all_labels, all_preds, normalize='true')
    classes = list(range(len(cm)))

    # Crear un DataFrame de la matriz de confusión
    df = pd.DataFrame(cm, index=classes, columns=classes)

    # Crear un heatmap de la matriz de confusión
    g = sns.heatmap(df, annot=True, cmap="Blues")
    g.set_yticklabels(g.get_yticklabels(), rotation=0)

    plt.title('Confusion matrix \n')
    plt.xlabel('Predicted label')
    plt.ylabel('True label')
    plt.autoscale()
    plt.show()
    print(f"test_loss_avg: {round(test_loss_avg,6)}\t test_accuracy: {round(test_accuracy*100,3)}%.")
    tf = time.time()
    h, m, s = time_elapsed(ti, tf)
    return test_loss_avg, test_accuracy

loss_test, acc_test = evaluar_red(dict_red1)
test_loss_avg: 0.135142	 test_accuracy: 97.64%.
Time Elapsed: 0.27 segs

En base a esto se puede ver que los resultados son excelentes, tanto para la clase ham (0) como para spam (1). Un punto importante a considerar es que las clases están desbalanceadas, por lo que al haber menos ejempos para spam, es más difícil aprender sobre esta. Quizás si se tuvieran más entradas de esta clase se podría obtener un mejor resultado con la misma arquitectura.

El accuracy en el conjunto de test es cercano al 98% y tiene buenos resultados para ambas clases, por lo que se tiene que la red sí logró aprender y generalizar. La gran cantidad de errores se deben a la clase spam y se estima que esto es por el desbalance de clases.

Algo se que se podría probar sería probar en el conjunto de entrenamiento a hacer DataAugmentation, haciendo padding cuando fuera necesario y cuando los textos fueran más grandes que el largo límite, se crearía un nuevo texto con lo que sobre. Se decidió no probar esto dado que los resultados son bastante buenos.

Función para tokenizar secuencias de texto:

In [28]:
def tokenize_review(test_review):
    test_review = test_review.lower() 
    test_text = ''.join([c for c in test_review if c not in punctuation])
    
    test_words = test_text.split()
    
    test_ints = []
    test_ints.append([vocab_to_int[word] for word in test_words])
    
    return test_ints

tokenize_review('hi we have an offer for u')
Out[28]:
[[104, 38, 17, 115, 374, 12, 6]]

Función para predecir si es spam o ham una secuencia de texto:

In [29]:
def predict(net, test_review, sequence_length):
      
    net.eval()
    
    test_ints = tokenize_review(test_review)
    # print(test_ints)
    
    seq_length = sequence_length
    features = padding_seq(test_ints[0], seq_length)
    
    feature_tensor = torch.from_numpy(features)
    
    batch_size = feature_tensor.size(0)
    
    h = net.init_hidden(batch_size)
    
    if(train_on_gpu):
      feature_tensor = feature_tensor.cuda()
      
    output, h = net(feature_tensor, h)
    
    pred = torch.round(output.squeeze())
    txt_print = 'Valor predicho: {:.4f}'.format(output.item())
    
    # print custom response based on whether test_review is pos/neg
    if(pred.item()==0):
      print(f'{txt_print}\t --->\t\033[94mHam\033[0m')
    else:
      print(f'{txt_print}\t --->\t\033[91mSpam\033[0m')

Ejemplo de uso:

In [30]:
predict(net,'hi we have an offer for u call back asap', 44)
Valor predicho: 0.0005	 --->	Ham

Función que realiza lo anterior con un ejemplo random del dataset de test:

In [31]:
def print_mod(texto, largo = 100):
    lineas = [texto[i:i+largo] for i in range(0, len(texto), largo)]
    print('\n'.join(lineas))
In [32]:
def probar_test(df, largo_secuencia = 44):
    idx = df.index.tolist()
    random_idx = random.choice(idx)
    random_row = df.loc[random_idx]

    texto = random_row['text']
    label = random_row['label']

    # Imprimir los valores seleccionados
    print('Texto:')
    print_mod(f'"{texto}"')
    if label == "ham":
        print(f'Label: \033[94m{label}\033[0m')
    else:
        print(f'Label: \033[91m{label}\033[0m')
    
    predict(net, texto, sequence_length = largo_secuencia)
    print(14*'---------')

Ejemplo de uso:

In [47]:
for i in range(5):
    probar_test(data_test, 44)
Texto:
"then u going ikea str aft dat"
Label: ham
Valor predicho: 0.0000	 --->	Ham
------------------------------------------------------------------------------------------------------------------------------
Texto:
"as a valued customer i am pleased to advise you that following recent review of your mob no you are
 awarded with a å£1500 bonus prize call 09066364589"
Label: spam
Valor predicho: 1.0000	 --->	Spam
------------------------------------------------------------------------------------------------------------------------------
Texto:
"want to funk up ur fone with a weekly new tone reply tones2u 2 this text wwwringtonescouk the origi
nal n best tones 3gbp network operator rates apply"
Label: spam
Valor predicho: 1.0000	 --->	Spam
------------------------------------------------------------------------------------------------------------------------------
Texto:
"you stayin out of trouble strangersaw dave the other day heåõs sorted nowstill with me bloke when u
 gona get a girl mrur mum still thinks we will get 2getha "
Label: ham
Valor predicho: 0.0000	 --->	Ham
------------------------------------------------------------------------------------------------------------------------------
Texto:
"88066 from 88066 lost 3pound help"
Label: spam
Valor predicho: 0.9990	 --->	Spam
------------------------------------------------------------------------------------------------------------------------------
In [34]:
# !jupyter nbconvert --to html CC6204_Tarea_2.ipynb

Conclusiones Generales¶

En general, se estima que se cumplieron con los objetivos de la tarea, logrando detectar spam y ham. Se comprobó la utilidad de las RNN vistas en clases y de las LSTM que aportan una mayor memoria a la red. Se evidencia que una de las partes más importantes es la limpieza del dataset, dado que casi siempre tienen datos erróneos o faltantes. Una de las partes más relevantes fue decidir el largo del padding, dado que esto puede mejorar o empeorar el desempeño del clasificador. Se utilizó un enfoque más estadístico para saber sobre el largo de los datos y tomar una decisión más informada. Por otro lado, se experimenta por primera vez con un embedding y se comprueba su utilidad y eficiencia en problemas de clasificación de texto.

Por último, se evalúa la eficacia de las redes neuronales recurrentes (RNN) en la clasificación de texto, destacando su capacidad para simplificar el proceso en comparación con las alternativas de modelos probabilísticos basados en gramática de NLP.